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Abstract 

Collisions between 248 Cm and 48 Ca are systematically investigated by time- 
dependent density functional calculations with evaporation prescription. De- 
pending on the incident energy and impact parameter, fusion, deep-inelastic 
and quasi-fission events are expected to appear. In this paper, possible fission 
dynamics of compound nuclei is presented. 

1 Introduction 

The synthesis of superheavy chemical elements 030 in the laboratory is one of the biggest challenges in 
nuclear physics. It is an attempt for clarifying the existence limits of all the chemical elements, as well 
as the completion attempt of the periodic table of chemical elements. We are concerned with heavy-ion 
collisions 

248 Cm + 48 Ca 

with different impact parameters in this paper. Let A and Z be the mass number and the proton number, 
respectively. The neutron number N is defined by N — A — Z, so that N/Z of 248 Cm and 48 Ca are 1.58 
and 1.40, respectively. If fusion appears, 296 Lv (= 296 116) with N/Z = 1.55 is produced. 

Fast charge equilibration (3) is expected to appear in low-energy heavy-ion reactions with an 
incident energy of a few MeV per nucleon. It provides a very strict limitation for the synthesis of 
superheavy elements. Actually, the N/Z of final product is not above nor below the N/Z of the projectile 
and the target (1.40 < N/Z < 1.58 in this case) in the case of charge equilibration, so that the proton- 
richness of the final product follows. Although the actual value of N/Z depends on the two colliding 
ions, its value for the merged nucleus tends to be rather proton-rich for a given proton number of the 
merged system. This feature is qualitatively understood by the discrepancy between the /3-stability line 
and the N = Z-line for heavier cases. Superheavy compound nuclei are very fragile and fission is a very 
frequent channel which leads to disintegration of the compound nuclei even at low excitation energies. 
In this paper, following the evaporation prescription shown in Ref. J4), a possible fission dynamics of 
compound nuclei is simulated based on the time-dependent density functional theory (TDDFT). 

2 Methods 

2.1 Treatment of the thermal property 

Self-consistent time-dependent density functional calculations are employed in this paper. TDDFT re- 
produces the quantum transportation due to the collective dynamics. In this sense, what is calculated by 
the TDDFT can be regarded as products after several 10 -21 s, which corresponds to a typical time-scale 
of low-energy heavy-ion reactions (1000 fm/c), as well as to the inclusive time interval of any collective 
oscillations such as giant dipole resonance, giant quadrupole resonance and so on. Meanwhile, ther- 
mal properties such as the thermal instability are not directly taken into account in TDDFT. Indeed, the 
Skyrme type interaction used in TDDFT (for example, see Ref. 0) is determined only from several 
densities. It is important that the most effective cooling effect arises from the emission of particles, and 
therefore it is expected that the break-up or fission of fragments including rather high internal excita- 
tion energy is suppressed in the TDDFT final products. The additional thermal effects leading to the 
break-ups of fragments should be introduced. 

Here is a fact that simplifies the treatment of thermal effects, that is, the difference of the time- 
scales. Different from the typical time scale of low-energy heavy-ion reactions, the typical time-scale 
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Fig. 1: (Colour online) Diagram of different reaction channels obtained in 48 Ca -f 248 Cm collisions by TDDFT 
calculations. The preferred reaction channels for different beam energies and impact paramters are given. The 
beam energies are all located above the Coulomb barrier which is 209.0 MeV. The results show fusion, deep- 
inelastic and elastic events. The difference between fusion-fission and quasi-fission is defined in this paper by 
whether the fission products satisfy "1.40 < N/Z < 1.58" or not. 



of the thermal effects is estimated by the typical time interval of collision-fission (fission appearing in 
heavy-ion collisions): several 10 -19 s. It is reasonable to introduce an evaporation prescription simply 
to the TDDFT final products. In this context the TDDFT final fragments have the meaning of products 
just after the early stage of heavy-ion reactions (several 10 -21 s). 

2.2 Evaporation prescription 

In complete fusion reactions the cross-section for the formation of a certain evaporation residue is usually 
given by three factors |[6l : 

crER(E cm ) = °Cp{E cm , J) x P C N(E cm , J) x Psv(E cm , J) (1) 
J 

where acp, Pcn and Psv mean the capture cross-section, the probability for the compound nucleus 
formation, and the probability for survival of the compound nucleus against fission. All three factors are 
functions of the centre-of-mass energy E cm and the total angular momentum J, where J can be related 
with the impact parameter. For light systems Pcn and Psv are about unity and ger ~ J2j g cp- But 
in superheavy systems the strong Coulomb repulsion and large angular momenta lead to small values 
of Pcn and Psv and therefore to the small cross-sections of the evaporation residues observed in the 
experiments. This means, different from light systems, it is necessary to introduce additional thermal 
effects for the superheavy element synthesis. First, acp is sufficiently considered in the TDDFT if we 
restrict ourselves to a sufficiently high energy exceeding the Coulomb barrier (cf. sub-barrier effects such 
as tunnelling are not taken into account in the TDDFT). Second, Pcn is fully considered in the TDDFT, 
which is a kind of mass equilibration also related to charge equilibration. Third, Psv whose relative 
time-scale is by no means equal to the former two probabilities is not satisfactorily considered in the 
TDDFT. This probability is much more related to thermal effects. Consequently, further consideration is 
necessary only for Psv as far as the energy above the Coulomb barrier is concerned. 

Several factors are included in Psv such as probabilities for fission of the compound nucleus, 
neutron-evaporation, proton-evaporation, deuteron-evaporation, alpha-particle-evaporation and so on. 
Probabilities of neutron and a-particle emissions are considered by 

PsV '•— (1 — Pn,evap)d ~ Pa,evap)- 

For the details of evaporation prescription, see our preceding research summarized in Ref. H. 

2.3 Fission dynamics 

For given energy and impact parameter, long-lived compound nucleus with certain excitation energies 
(compared to the ground state) are obtained. Fission dynamics is obtained by additional TDDFT cal- 
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Fig. 2: (Colour online) Time evolution of 248 Cm -f 48 Ca with a fixed incident energy 268 MeV in the centre- 
of-mass frame, where a box is fixed to 48x48x24 fm 3 . Snapshots at 2.3xl0~ 22 s, 8.7xl0~ 22 s, 15.1 xlO -22 s 
and 21.5 x 10 -22 s are shown, where a snapshot at 21.5 x 10 -22 s is not shown only for b = 10 fm. Fusion, deep- 
inelastic and elastic events appear depending on b. In particular, the life-time of composite nucleus for b = 6 and 
8 fm is the order of x 10 -21 s. 
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culations, where no initial velocity is given. Let us take a residual nucleus of mass number Ar, proton 
number Zr and excitation energy Er. Consider the binary fission: 

Ar Zr z x + M Z 2 , (2) 

where Ar = Y^h=i A% and Zr = Y^=i First, choose A{ and Z{ to determine the heavy-ion reaction 
being considered. A configuration of the two nuclei at a distance Rq is prepared as an initial state 
of additional calculation. Second, choose Rq such that the excitation energy agrees with that of the 
compound state found in the collision (Fig. 0), where the energy can be different depending on Rq. Note 
here that the TDDFT is a theory in which the total energy is strictly conserved, so that the total energy 
is conserved during the presented fission process. Third, the initial many-body wave function, which 
is given as a single Slater determinant, consists of single wave functions of two different initial nuclei, 
where a set of single wave functions are orthogonalized before starting TDDFT calculations (cf. the 
Gram-Schmidt orthogonalization method). In this way a configuration at the same excitation energy but 
closer to fission can be obtained; for given Ar, Zr and Er. The distance Ro is uniquely determined for 
fixed M Z X and M Z 2 . 



3 TDDFT calculations for fission dynamics 

Time-dependent density functional calculations with a Skyrme interaction (SLy6 [7]) are carried out in a 
spatial box of 48 x 48 x 24 fm 3 with periodic boundary condition. The unit spatial spacing and the unit 
time spacing are fixed to 1.0 fm and 2/3 x 10 -24 s, respectively. 

The initial positions of 248 Cm and 48 Ca are fixed to (0, 6, 0) and (-15,0,0), respectively. The initial 
248 Cm is almost spherical the diameter for x, y, and z directions are 19 fm, 19 fm and 18 fm, respectively. 
248 Cm (the right hand side on x — y -plane) does not have initial velocity on the frame, while the initial 
velocity parallel to the x-axis is given to 48 Ca (the left hand side on x — ?/-plane). The systematic 
results of TDDFT calculations for a given incident energy (268 MeV) are summarized in Fig. [TJ These 
results, which include many fusion events, provide a quite optimistic view for producing superheavy 
elements. However, in comparison with experiments, the corresponding fusion cross-section of those 
low-energy heavy-ion reactions is too high to believe. Consequently, although these TDDFT results are 
still legitimate to show products just after the early stage of heavy-ion reactions, it is necessary to take 
into account Psv in order to have comparable results to experiments. 

In case of the incident energy 268 MeV, the pure TDDFT results show the following reactions 
(Fig.©: 

248 Cm + 48 Ca ^ 296 Lv 
248 Cm + 48 Ca ^ 247 Cf + 49 Ar 
248 Cm + 48 Ca ^ 246 Bk + 50 R 
248 Cm + 48 Ca 248 Cm + 48 Ca 

for b = 4, 6, 8, 10 fm, respectively. If we take into account the neutron and alpha emissions, they become 

Ds + 3a + 4n 
Pu + 47 Ar + 2a + 3n 
Np + 48 K + 2a + 3n 
Cm + 48 Ca 

for b = 4, 6, 8, 10 fm, respectively. 

Concerning the fission dynamics, here we take the case of b = 4 fm. For instance we consider the 
symmetric fission for 296 Lv, which corresponds to the pure TDDFT product. The distance Rq =11.8 fm 
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Fig. 3: (Colour online) Symmetric fission dynamics from the compound nucleus, which corresponds to the 
heavy-ion reaction shown in Eq. ©, is simulated. A calculation box is fixed to 48x48x24 fm 3 . Snapshots at 
0.0 x 10" 22 s, 8.3 x 10" 22 s, 16.6 x 10" 22 s and 24.9 x 10~ 22 s are shown. The distance R , which is deduced from 
the excitation energy, is equal to 11.8 fm. Note that the orthogonalization is applied to the initial state before 
starting TDDFT calculations, because the two initial nuclei are slightly overlapped. 



is deduced from the excitation energy of the compound nucleus ( Lv), where trial TDDFT calculations 
with several Ro were performed to identify Ro. The fission dynamics 

248 Cm + 48 Ca 296 Lv 148 Ce + 148 Ce (3) 

is shown in Fig. [5J The magnitude of Ro is related to the difficulty of fission, as well as the total time- 
scale of fission. In addition Rq becomes larger for lower incident energies. Note that Ro can be larger 
than the touching distance in which case the fission is suggested to be impossible. In this sense this 
method is applicable to the fission appearing in heavy-ion collisions. Although the duration of fission 
shown in Fig. [3] is quite short (similar to the typical duration time of low-energy heavy-ion reactions), 
the additional time is necessary to realize the initial state. The total duration time is expected to be 
significantly longer than 10 -21 s, because the initial state shown in Fig. [51 which cannot be realized in 
the standard TDDFT at the least, cannot be easily realized. 

4 Summary 

A procedure of obtaining fission dynamics of excited compound nuclei has been presented, where the 
reproduction of microscopic fission dynamics was a long standing problem in nuclear theory (among 
a few preceding works on microscopic fission dynamics, see J. W. Negele O). The presented method 
allows us to have a self-consistent treatment of fission dynamics. The obtained dynamics treating the 
excited states under the strict total energy conservation is diabatic, which is essentially different from the 
adiabatic fission dynamics. 

For the proposed method it is remarkable that the initial condition of fission dynamics is uniquely 
determined without having any intentional settings; i.e., it is automatically determined for a given set of 
Ar, Zr, Al Zi, A2 Z2 and Er. As is seen in the comparison between Figs. [2] and [51 the presented fission 
dynamics is hidden if only pure TDDFT calculations are utilized. In this sense a new point of this method 
is to choose an ideal configuration (which is rarely realized in most cases if the excitation energy is not so 
high) as the initial state. As a matter of cause, a different choice of Al Z\ and A<1 Z2 brings about different 
fission dynamics even from an identical compound nucleus. It is also worth noting that this method with 
taking into account many different initial conditions can be utilized to distinguish whether fission arising 
from the collective dynamics can appear or not. Even though the proposed procedure might not be the 
ultimate solution for investigating fission dynamics, the obtained dynamics actually extracts important 
aspects of fission dynamics (for example, see the time evolution of the neck). 

It was suggested by the calculation shown in Fig. [5] that symmetric fission of the compound nu- 
cleus ( 296 Lv) is possible in the collision: 248 Cm + 48 Ca with an incident energy E =268 MeV. The 
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fission process itself takes only a few 10 s, but the preformation of the initial state might take signifi- 
cantly longer. The estimation of the total duration time for fission process is a future problem. 

This work was supported by the Helmholtz alliance HA216/EMMI. The authors thank Prof. J. A. 
Maruhn for reading this manuscript carefully. 
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